exp(-max(0, x - 31)^2 / (2 * -((365*3)^2 / (2*ln(0.5)))))* log10(y + 2) plot 0 <=x <= 1000, 1 <= y <= 1000